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Abstract 


A calculation* al model is derived for use in estimating Solar cosmic 
ray exposure to critical body organs in low Earth orbit at the center of 
a large spherical shield of fixed thickness . The effects of the Earth’s 
geomagnetic field and the astronauts’ self-shielding are evaluated ex- 
plicitly . The geomagnetic field model is an approximate tilted eccentric 
dipole with geomagnetic storms represented as an impressed uniform 
field . The storm field is related to the planetary geomagnetic index K p . 
The Shuttle geometry is introduced into the resultant computer code us- 
ing the Shuttle mass distribution surrounding two locations on the flight 
deck . The Shuttle is treated as pure aluminum and the astronaut as soft 
tissue . Short-term average fluence over a single orbit is calculated as a 
function of the location of the lines of nodes or long-term averages over 
all lines of nodes for a fixed inclination. 


Introduction 

Solar cosmic rays observed in low Earth orbit first 
passed through the Earth’s magnetic field. Those 
particles that are able to penetrate the geomagnetic 
field must further penetrate the walls of the space- 
craft before exposing the human occupants* As a 
result of interactions in the vehicle structure and the 
bulk tissues of the astronauts 1 bodies, the composi- 
tion of the rays is greatly altered. Any reasonable 
estimate must account for geomagnetic effects, the 
atomic and nuclear interactions, and the spacecraft 
and human body geometry. 

During times of increased solar activity, small 
amounts of the solar plasma are ejected into inter- 
planetary space. When this plasma interacts with the 
Earth’s magnetic field, large distortions of this field 
result in geomagnetic storms. Since the energetic 
solar flare particles often arrive during such geomag- 
netic disturbances, the penetration of the energetic 
protons into the magnetospheric cavity can be vastly 
different from that seen under quiet conditions. 

In a previous report (ref. 1), we derived a simple 
model for estimating exposure to the critical organs 
of an astronaut within a large spherical shield using a 
centered-tilted dipole model of the geomagnetic field 
and an approximate geomagnetic storm field as an 
impressed uniform field. In the present work, ap- 
proximate expressions are derived for a tilted eccen- 
tric dipole model. Geomagnetic model parameters 
are used to better approximate the recently calcu- 
lated worldwide vertical cutoffs of Shea and Smart 
(ref. 2) for the 1980 magnetic field parameters than 
in the previous study (ref. 1). The average transmis- 
sion factors of the approximate model agree well with 
the transmission factors of Smart and Shea (ref. 3). 
The approximate orbit-averaged transmission factors 


are found as a function of location of the ascending 
line of nodes, and dose is given for the worst-case 
exposed orbit for various geomagnetic storm condi- 
tions. A comparison is made with earlier work. 

A necessary requirement for space mission analy- 
sis is a reliable method of calculating the anticipated 
dose distribution in the human body. In the case 
of extraterrestrial radiations, a model is required for 
evaluating the effects of the geomagnetic field in ad- 
dition to the interaction with the spacecraft struc- 
ture and the human body. Described herein are ad- 
ditional computer subroutines which calculate the 
proton dose within the Shuttle averaged over five 
major segments (upper limbs, lower limbs, upper 
trunk, lower trunk, and skull) of the blood-forming 
organs (BFO), ocular lens, and skin by treating the 
human body and vehicle geometry in detail (refs. 4 
and 5) but assuming isotropicity of the incident pri- 
mary particles. The calculation uses an approximate 
form of transport theory in which nuclear star ef- 
fects are incorporated (refs. 1 and 6). The output 
of the program is in terms of physical dose and dose 
equivalent. 

Geomagnetic Effects on Orbital 
Environment 

Charged particles arriving at some location within 
the geomagnetosphere are deflected by the Lorentz 
force evxB which prevents penetration for some 
directions of incidence and some energies. Such 
phenomena were extensively studied by Stormer 
(ref. 7) for a dipole magnetic field which provides 
the basis for classifying the orbital trajectories of 
charged particles arriving at some location within 
the field. As a part of Stormer ’s theory, there ex- 
ist allowed trajectories with no connection to asymp- 
totic trajectories now recognized as trapping regions 



associated with Van Allen radiation. Numerical solu- 
tions to the charged particle equations of motion in 
more realistic geomagnetic field models were intro- 
duced by McCracken (ref. 8) and further advanced 
by Smart and Shea (refs. 2, 3, and 9). Our pur- 
pose here will not be to supplant the vastly detailed 
numerical work but to seek a simple analytic form to 
reasonably approximate the more general numerical 
solutions. 

The geomagnetic field can be reasonably 
approximated by a tilted dipole with moment 
M = r 3 31 500 nT displaced from the Earth’s center 
by 430 km or 0.0G8r e where r c = 6378 km. The 
tilt angle is 11.7° at 69° W longitude. The magnetic 
quadrupole contributions are then about 10 percent 
at the surface and decrease to 5 percent at 2r e . 
Higher-order moments are even smaller. The mo- 
tion of charged particles in the geomagnetic field was 
studied extensively by Stormer. We outline his meth- 
ods herein. In spherical coordinates, Stormer showed 
that the azimuth (p is an ignorable coordinate pos- 
sessing an integral for the particles trajectory such 
that 
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where m is the mass of the particle, Ze is the charge. 
v is the velocity, c is the velocity of lighte r is radial 
distance from the center of the field, 9 is colatitude, 7 
is an integration constant, and uj is the angle between 
the velocity vector and the azimuthal direction. The 
allowed Stormer regions consist of the space for which 
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Further analysis of the condition in equation (2) 
shows stable trapping regions as well as the 
Stormer main cone of transmission given for 
7 = 2mv(ZeM/mvc)^^. The Stormer main cone is 
given (ref. 10) by the solid angle element 


4n reaching its half-value at uj — tt/2 including an- 
gles up to the vertical direction. The vertical cutoff 
model is expressed as 


fi « 4^ U [R - R c (X m )] (5) 

where the vertical cutoff rigidity from equation (4) is 

M 4 

R C Om) = -^2 COS Am ( 6 ) 


and U(x) is the unit step function. 

Not included in the above formalism are those 
trajectories which are cut off by the shadow cast 
by the solid Earth. The fraction of the solid angle 
covered by the Earth’s shadow is estimated, assuming 
that the curvature of the local trajectories is large 
compared with the Earth’s radius (ref. 10). Then 
the solid angle fraction is 
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The corrected solid angle for the vertical cutoff model 
is then 

fi = fi s h U [i? - Rc (A m )] (8) 

which leaves the local solid angle open to transmis- 
sion of charged particles of rigidity R at altitude r 
and geomagnetic latitude A m . 

Spacecraft in low Earth orbit (LEO) are typically 
in circular orbits which simplifies the analysis. The 
orbit plane is inclined with respect to the equato- 
rial plane. Since the Earth’s angular momentum 
(spin) and the spacecraft orbital angular momentum 
are conserved, the angle between them is fixed and 
equal to the inclination angle i. The magnetic axis 
rotates with the Earth and therefore precesses about 
the rotational axis within a 24-hour period. The geo- 
graphic location of the ascending node likewise moves 
around the geographic equator every 24 hours. The 
inclination of the orbit plane relative to the magnetic 
axis i m likewise is periodic. If 77 is the geographic co- 
ordinate of the ascending node line, then 


fi = 27r(l + cos uj ) (3) 

which contains the allowed directions of arrival for 
particles of rigidity R (momentum per unit charge) 
given by 


cos i m — cos i cos 9 m -j- sin i sin 9 m cos (77 - (p m - 90°) 

( 9 ) 

where 9 rn and <p m are the magnetic north pole colati- 
tude and longitude. The average transmission factor 
around this orbit F is then 


R = 


M 
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Henceforth we replace the colatitude 9 by the mag- 
netic latitude X m and note that fi varies from 0 to 


F(R,i. n)=%± f 7 ” U{R-R c (X)]dX= i ^ULn sh (10) 

l m J 0 im 

where X m is the magnetic latitude with cutoff at R 
as given by 

R C (\m)=R (11) 
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We note that i m goes through a maximum and 
minimum orbit corresponding to rj m ax = ~ 90 

and 77 m i n = (pm + 90° for which i m = i + 0 m and 
i m = \i - 9 m |, respectively, as we have shown else- 
where (ref. 1). 

We may also calculate the long-term average over 
many days of orbits by averaging equation (10) over 
the node angle 77 as 
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where 


cos i m = cos ^ cos + sin i sin 6 m cos (p (13) 


and 
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Equation (12) may be rewritten as 
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where the last integral is approximated by a numer- 
ical quadrature. The results of equation (15) are 
compared with the numerical calculations of Smart 
and Shea for 400-km (216-n.mi.) orbits at several in- 
clinations in figure 1 for this centered dipole field 
model with 9 m = 11.7° (tilt angle) and longitude 
4hn = -69° (69° W). 

An important correction to the centered dipole 
field is the displacement of the geomagnetic dipole 
430 km (232 n.mi.) from the Earth’s center. Unfor- 
tunately, the formalism is very complicated, since the 
distance r from the dipole center is no longer constant 
even for a circular orbit. The offset dipole decreases 
the cutoffs in the Atlantic hemisphere defined by the 
meridinal plane normal to the tilt direction and in- 
creases the cutoffs in the remaining hemisphere over 
the Pacific. We define two cutoff functions for cen- 
tered dipole fields as 
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(c) Orbit inclination, 40°. (d) Orbit inclination, 50°. 


Figure 1. Dipole maximum and minimum cutoff model and 
numerical simulation of exact geomagnetic field model. 


where j = A, 8 a = 593 km (320 n.mi.) for the At- 
lantic hemisphere; j — P, 8p = —504 km (-272 n.mi.) 
for the Pacific; A m is the usual magnetic latitude 
which depends on the hemisphere; and r is the geo- 
centric radius of the orbit. The value 14.9 GV 
is found from the value of the dipole moment of 
r 3 31 500 nT, and 8j were chosen to match the mini- 
mum equatorial cutoff in the Atlantic region of Shea 
and Smart (ref. 9) and the maximum cutoff in the 
Pacific. The calculation of the orbit average trans- 
mission factor is as before except that the two hemi- 
spheres are considered separately as 


F(R,i,r]) — 2^ 2 sh 
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im 
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where 

Ra(W = R ( 18 ) 

and 

Rp(\p) = R ( 19 ) 

Similarly, the long-term average of equation (12) is 
extended to each hemisphere as 



where the integrals are evaluated as described for 
equation (15). The average transmission factors 
of equation (20) are compared with the numerical 
calculations of Smart and Shea (ref. 3) in figure 2. 
The pole’s tilt angles (ref. 11) are given in table 1 
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along with suitably chosen longitudes and are shown 
in relation to the vertical cutoff rigidities of Smart 
and Shea in figures 3 and 4. 



(a) Orbit inclination, 40°. 



(b) Orbit inclination, 50°. 

Figure 2. Offset dipole model average transmission factors 
with detailed calculations of Smart and Shea (ref. 3) and 
maximum and minimum transmission factors. 

Table 1. Geographic Locations of Offset Poles in 
Present Calculations 


Magnetic 

Longitude, 

Tilt, 

pole 

0m, beg 

0m i deg 

North 

-G9 

16 

South 

121 

22 


During times of intense solar activity, the solar 
plasma emitted in solar flares and subflares advances 
outward and arrives at 1 AU from the Sun. If the 


Earth is locally present, the plasma interacts with 
the geomagnetic field in which the plasma pressure 
performs work on the local geomagnetic field. The 
initial impact produces hydromagnetic waves caus- 
ing a general increase in geomagnetic intensity. As 
plasma flow is established, it generates large electric 
ring currents and a corresponding impressed mag- 
netic storm field. In the initial phase (hydromagnetic 
wave), the storm field is parallel to the equatorial 
field after which the storm field reverses in the main 
phase of the storm caused by ring currents within the 
magnetopause and opposes the quiet field, to cause a 
net decrease of the field strength. The main phase is 
followed by slow recovery to the quiet field conditions 
(ref. 11). 

The magnetic storm model used here assumes a 
uniform magnetic field impressed on the normal quiet 
field (ref. 10). The storm field strength can be found 
from the change in the horizontal field component 
around the geomagnetic equator. We represent this 
field by H s t- Typical values of H s t in the main phase 
range from substorm values -10 nT to severe storms 
with —500 nT. On rare occasions, for very intense 
storms, the storm field exceeds -1000 nT. 

Magnetic disturbances have been observed for 
many years, and various classification schemes for 
such disturbances have been proposed. The plane- 
tary magnetic index K p is based on magnetometer 
measurements of 12 stations worldwide. The K p in- 
dex is related to a derived planetary index a p and 
storm field strength by Bartels (ref. 11) given in 
table 2. 

Table 2. Relation of Magnetic Indices to Magnetic 
Storm-Field Strength 


Kp 

dp 

IflstUT 

0 

0 

0 

1 

4 

8 

2 

7 

14 

3 

15 

30 

4 

27 

54 

5 

48 

96 * 

6 

80 

160 

7 

132 

264 

8 

207 

414 

9 

400 

800 


The vertical cutoff rigidity as given by equa- 
tion (16) is further modified to approximate the ef- 
fects of geomagnetic disturbances. It was shown by 
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Kuhn, Schwamb, and Payne (ref. 10) that the appro- 
priate equation is 
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for the centered dipole field. In the context of our 
approximation of the offset-tilted dipole field, we get 
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This vertical cutoff replaces equation (16) and applies 
to storm conditions. Note that the cutoff is zero 
whenever the result of equation (2 2 ) is negative. The 
corresponding transmission factor on the worst case 
orbit (77 £3 211°) is shown in relation to the quiet 
field average transmission factors of Smart and Shea 
(ref. 3) in figure 5. 


Dose Estimation 

In passing through tissue, energetic protons inter- 
act mostly through ionization of atomic constituents 
by the transfer of small amounts of momentum to 
orbital electrons. Although the nuclear reactions are 
far less numerous, their effects are magnified because 
of the large momentum transferred to the nuclear 
particles and the struck nucleus itself. Unlike the 
secondary electrons formed through atomic ioniza- 
tion by interaction with the primary protons, the ra- 
diations resulting from nuclear reactions are mostly 
heavy ionizing and generally have large biological ef- 
fectiveness. Many of the secondary particles of nu- 
clear reactions are sufficiently energetic to promote 
similar nuclear reactions and thus cause a buildup of 
secondary radiations. The description of such pro- 
cesses requires solution of the transport equation. 
The approximate solutions for the transition of pro- 
tons in 30- cm- thick slabs of soft tissue for fixed inci- 
dent energies have been made (ref. 6). The results of 
such calculations are dose conversion factors for re- 
lating the primary monoenergetic proton fluence to 
dose or dose equivalent as a function of position in a 
tissue slab. 

Whenever the radiation is spatially uniform, the 
dose at any point x in a convex object may be 
calculated (ref. 12) by 



Rigidity, GV 
(a) Orbit inclination, 30°. 



Rigidity, GV 
(b) Orbit inclination, 40°. 



Rigidity, GV 


(c) Orbit inclination, 50°. 


D{x) = r f R u [z x (n), E]4>{Q, E) dn dE (23) 

J 0 J SI 


Figure 5. Maximum transmission factor for various storm 
fields and quiet time average transmission of Smart and 
Shea. 7] ~ 211°. 


where R n {z, E) is the dose at depth 2 for normal in- is the local differential proton fluence along direction 

cident protons of energy E on a tissue slab, <p{Q, E) 0, and z x (ft) is the distance from the boundary 
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along ft to point x . It has been shown that equa- 
tion (23) always overestimates the dose but is an ac- 
curate estimate when the ratio of the proton beam 
divergence due to nuclear reaction to the bodies’ ra- 
dius of curvature is small (ref. 12). Equation (23) is 
a practical prescription for introducing nuclear reac- 
tion effects into calculations of dose in geometrically 
complex objects such as the human body. The main 
requirement is that the dose conversion factors for a 
tissue slab be adequately known for a broad range of 
energies and depths. The dose conversion factors for 
tissue were derived in reference 6, and a correction 
for an aluminum shield was found in reference 1. The 
spacecraft geometry is taken as an aluminum sphere 
of large radius. 


Method for Shuttle Geometry 

In the previous section, the calculation of astro- 
naut exposure in the center of a large aluminum 
sphere of arbitrary thickness was derived for a specific 
orbit with either the quiet geomagnetic field or with 
a geomagnetic disturbance. We denote that result by 
D sp h {t s )j and it has a different value for each criti- 
cal organ for which exposure within the aluminum 
sphere of thickness t s is evaluated. W ithin the con- 
text of assumed isotropic radiation, the exposure at 
some location within the Shuttle is 

rCfO 

D = / D sph (t s ) f(t 8 ) dt s (24) 

where f{t s ) describes the mass distribution of the 
Shuttle structure assumed to be aluminum about 
that particular location. Physically, f(t s ) dt s is the 
solid angle fraction for which the areal density to the 
Shuttle surface lies between t a and t s + dt s • The 
cumulative distribution of areal density is given by 


and 


fiits) 


0.303 

IT 

0.147 

ts 


(1 <ts< 6) 

(6 < t a < 132) 


(27) 


where the functions are understood to be zero outside 
of specified ranges. Formulas (24), (26), and (27) are 
used in conjunction with the methods described in 
previous sections to estimate Shuttle exposure in the 
two locations referred to in reference 5 as dosimeter 
location numbers 1 and 2 which comprise the most 
and the least shielded locations in the Shuttle crew 
compartment, respectively. The method can be eas- 
ily expanded to include more astronaut organs and 
other Shuttle locations. 



Figure 6. Mass distribution of two locations on Shuttle flight 
deck. 


Ffj [t s ) = f f{t s ) dt s (25) 

JO 

and is shown for two locations in the Shuttle (ref. 5) 
in figure 6. Also shown in figure 6 are the following 
approximate functions: 



0.176 

ts 

(1 < ts < 2) ' 

flits) = < 

0.113 

ts 

(2 <t s < 20) > 


0.353 
k t$ 

(20 < t s < 120) ^ 


Results 

The maximum exposure limits in force for the 
Space Station Freedom (ref. 13) are shown in ta- 
ble 3. The dose and dose equivalent to critical body 
organs for an aluminum shield 1 g/cm thick are 
shown in tables 4 through 7 for various storm con- 
ditions (Hst)- The exposures are shown for a worst- 
exposed orbit ( r) = 211°) and the average over all 
r]. The average is shown, since average transmission 
factors are calculated by several groups, and one may 
be tempted to use the transmission factor appropri- 
ate for the galactic cosmic ray background (refs. 14 
and 15). It is clear from the results in tables 4 
through 7 that such use of average cutoffs provides 
exposure estimates which could be too small by a 
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factor of 2 to 10. Such an underestimate is clearly 
unacceptable. Furthermore, comparing the current 
dose estimates with values for transmission factors 
derived for a tilted concentric dipole field (ref. 1), we 
see that the dose values of the current field model are 
a factor of 3 to 5 higher. The eccentric field had two 
effects which lower the cutoffs. The offset displaces 
the South magnetic pole to lower latitudes and low- 
ers the geomagnetic cutoff values over the Atlantic. 
The methods derived herein allow evaluating expo- 
sures as a function of the location of the line of nodes 
and should provide acceptable estimates of exposure. 

Table 3. Ionizing Radiation Exposure Limits for Space 
Station Freedom Astronauts 



Dose equivalent, cSv, for — 

Exposure 

interval 

Skin 

Eye 

Blood-Forming 

organs 

30 days 

150 

100 

25 

Annual 

300 

200 

50 

Career 

600 

400 

“100-400 


“Dependent on gender and age at initial exposure. 


Observing the levels of exposure in low inclina- 
tion orbits ( i « 30°), a significant exposure could 
clearly occur if particle arrivals coincided with a large 
magnitude (K p ~ 9) magnetic disturbance. On the 
basis of the present analysis, a more in-depth study 
of potential solar flare exposure of the Space Station 
Freedom seems warranted. Such a study should in- 
clude a review of the history of major geomagnetic 
disturbances in proximity to solar particle events, a 
review of alternate geomagnetic storm models, and 
a review of the specific Space Station Freedom shield 
geometry. 


The exposure for Shuttle flight in a 400-km 
(216-n.mi.) orbit with a 50° inclination are showm 
in table 8 for the February 23, 1956, solar event 
spectrum as compiled by Foelsche et al. (ref. 16). 
A magnetic storm was assumed to be in progress 
with an impressed field of —100 nT. The results 
shown in table 8 are for the long-term, average geo- 
magnetic cutoffs, since these are directly comparable 
with the work of other geomagnetic models (refs. 14 
and 15). We note, however, that actual exposure 
could be greatly different depending on the location 
of the line of nodes at the time of arrival of the high- 
energy flare particles. 

Concluding Remarks 

The present code should be useful for assessing 
the potential exposure in low Earth orbit missions if 
the Solar spectrum and state of the geomagnetic field 
arc known. - 

It was shown for low inclination orbits ( i & 30°) 
that a significant exposure could occur if the parti- 
cle’s arrival coincided with a large magnitude ( K p « 
9) magnetic disturbance. On the basis of the present 
analysis, a more in-depth study of potential solar 
flare exposure of Space Station Freedom seems war- 
ranted. Such a study should include a review of the 
history of major geomagnetic disturbances in prox- 
imity to solar particle events, a review of alternate 
geomagnetic storm models, and specific Space Sta- 
tion Freedom shield geometry. 

NASA Langley Research Center 
Hampton, VA 23665-5225 
April 17, 1990 



Table 4. Skin Dose Behind Aluminum Shield 1 g/cm 2 Thick During February 25, 1956, November 12 13, 1960, and August 4, 1972, 

Events With Various Storm Fields 



Table 5. Skin Dose Equivalent Behind Aluminum Shield 1 g/cm 2 Thick During February 25, 1956, November 12-13, 1960, and August 4, 1972 ; 

Events With Various Storm Fields 





Table 6. BFO Dose Behind Aluminum Shield 1 g/cm 2 Thick During February 25, 1956, November 12-13, 1960, and August 4, 1972, 

Events With Various Storm Fields 


Orbit 

inclination, 

deg 

BFO dose, cGy, during - 

Feb. 1956 for H s t, nT, of- 

Nov. 1960 for H s t, nT, of 


-100 

-500 

-900 

-100 

-500 

-900 

-100 

-500 

-900 

30, max 

< 0.1 

1.0 

3.3 

0 

.9 

4.4 

0 

1.8 


30, avg 

< .1 

.2 

.8 

0 

.1 



.2 

2.3 

40, max 

1.4 

4.9 

6.8 

1.3 

6.7 

9.5 

2.5 

16.0 


40, avg 

.3 

1.5 

2.5 

.2 

2.0 

3.4 

.4 

4.7 

8.8 

50, max 

4.6 

7.6 

9.3 

6.2 

11.0 

13.0 



33.0 

50, avg 

1.5 

3.2 

4.4 

1.8 

4.4 

6.2 

4.2 

11.0 

16.0 


Table 7. BFO Dose Equivalent Behind Aluminum Shield 1 g/cm 2 Thick During February 25, 1956, November 12 13, I960, and August 4, 1972, 

Events With Various Storm Fields 


Orbit 

inclination, 

deg 

BFO dose equivalent, cSv, during — 

Feb. 1956 for H st , nT, of 

Nov. 1960 for H s t, nT, of — 

Aug. 1972 for H st , nT, of — 

-100 

-500 

-900 

-100 

-500 

-900 

-100 

-500 

-900 

30, max 

< 0.1 

1.7 

5.1 

0 

1.3 

6.1 

0 

2.6 

15.0 

30, avg 

< .1 

.3 

1.2 

0 

.2 

1.3 

0 

.3 

3.3 

40, max 

2.3 

7.6 

10.0 

1.8 

9.2 

13.0 

3.6 

23.0 

34.0 

40, avg 

.5 

2.3 

3.9 

.3 

2.7 

4.7 

.5 

6.7 

12.0 

50, max 

7.3 

12.0 

14.0 

8.5 

15.0 

18.0 

21.0 

38.0 

48.0 

50, avg 

2.3 

5.0 

6.9 

2.5 

6.1 

8.6 

6.0 

15.0 ! 

22.0 


Table 8. Human Exposure at Two Locations in Shuttle Crew Compartment for February 23, 1956, 
Event With H s t = -100 nT at 400 km and Orbit Inclination of 50° 



Exposure in — 


BFO 

Skin 

Lens 

Location 

cGy 

cSv 

cGy 

cSv 

cGy 

cSv 

1 

2.6 

4.9 

4.0 


4.1 

7.6 

2 

3.4 

5.9 

6.0 


6.0 

10.0 


10 
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